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Understanding the mieroscopic processes afFeeting the bulk thermal eonduetivity is erucial to de- 
velop more efhcient thermoelectric materials. PbTe is currently one of the leading thermoelectric 
materials, largely thanks to its low thermal conductivity. However, the origin of this low thermal 
conductivity in a simple rocksalt structure has so far been elusive. Using a combination of inelastic 
neutron scattering measurements and first-principles computations of the phonons, we identify a 
strong anharmonic coupling between the ferroelectric transverse optic (TO) mode and the longitu- 
dinal acoustic (LA) modes in PbTe. This interaction extends over a large portion of reciprocal space, 
and directly affects the heat-carrying LA phonons. The LA-TO anharmonic coupling is likely to 
play a central role in explaining the low thermal conductivity of PbTe. The present results provide 
a microscopic picture of why many good thermoelectric materials are found near a lattice instability 
of the ferroelectric type. 



Thermoelectric materials are of intense interest for 
energy applications, because they can transform heat, 
otherwise lost to the environment, to produce poten- 
tially useful electricity [1-3]. The thermoelectric figurc- 
of-merit, zT = S'^aT/n, which determines the maximum 
efBciency of the conversion process, is determined by the 
Seebeck coefficient, S, the electrical conductivity, a, and 
the thermal conductivity, k [3]. This formula expresses 
the fact that the electrical energy produced is in the 
form of a current driven by the thermoelectric voltage, 
which, under open circuit conditions, is given by —SAT, 
while heat conduction and electrical resistance are para- 
sitic. Thus, efficient thermoelectric materials must have 
a low thermal conductivity, in order to preserve the tem- 
perature gradient exploited to generate a voltage. The 
thermal conductivity contains electrical (kci) and lattice 
(K^iat) components. The electrical part is related to the 
electrical conductivity and therefore reducing Kjat is cen- 
tral to obtaining high performance thermoelectrics. A 
powerful approach consists in designing materials with 
defects and micro/nanostructural features that hinder 
the propagation of phonons (the atomic vibrations that 
carry most of the heat in semiconductors), while pre- 
serving a [2, 4]. However, the difficulty of fabrication 
of materials with such tailored nanostructures, and their 
thermodynamic stability at high-temperature, can be a 
concern for real-world applications. Thus, thermoelectric 
materials with an intrinsically low Kiat in bulk form are 
of very high interest. 

PbTe is one of the leading thermoelectric materials in 

the temperature range 400 to 800 K [5] , with a very low 
Kiat, as well as a large Seebeck coefficient when appro- 
priately doped, and a good electrical conductivity [1, 6]. 
However, the origin of the low Kiat is not well under- 
stood. PbTe crystallizes in the rocksalt structure, a sim- 
ple high-symmetry structure a priori not associated with 
low thermal conductivity. The bulk K\at in PbTe is indeed 
surprisingly low, with Kiat = 2Wm~^K~-'^ at 300 K in 
single crystalline samples [7], and similar values in poly- 
crystalline samples [8]. Thus, a perfect single crystal of 



PbTe is about five times more resistant to heat conduc- 
tion than Sio.yGeo.a (k = 10Wm~^K~^), a thermoelec- 
tric alloy with strong random mass disorder. Generally, 
Kiat is expected to be low in disordered alloys, materials 
with nanoscale inhomogeneities, or crystals with com- 
plex unit cells containing loosely bound atoms "rattling" 
in cages [3], all features that can lower k by scattering 
phonons [9-12]. Clearly, another mechanism, imknown 
until now, must be responsible for the low k in rocksalt 
PbTe. 

Thanks to recent developments in instrumentation, in- 
elastic neutron scattering (INS) instruments, in particu- 
lar time-of-flight spectrometers at spallation sources, are 
now able to measure the entire four-dimensional scatter- 
ing function, S{Q,E), which contains the full informa- 
tion about the microscopic dynamics of materials. These 
measurements provide much more detail than was previ- 
ously obtainable with usual triple-axis neutron scatter- 
ing experiments [13, 14]. In parallel, powerful ab initio 
computational techniques can be used to calculate the 
phonon dispersions, and can be compared directly with 
experiments. 

Here, we present a detailed mapping of the four- 
dimensional ^(Q,^) of PbTe based on inelastic neu- 
tron scattering measurements on single-crystalline PbTe, 
and results of ab initio phonon calculations using den- 
sity functional theory (DFT). While our measurements 
are in general agreement with previous reports [13-15], 
our more detailed investigations also reveal several new 
key features that directly relate to the origin of the low 
K. In particular, our measurements of the tempera- 
ture and wavevector dependences of the TO mode reveal 
previously unreported anomalies. In particular, we ob- 
serve the signature of a very strong and extended anhar- 
monic LA-TO coupling. The LA-TO coupling leads to 
an avoided-crossing behavior in the dispersions, as well as 
an anomalous lowering and damping of the LA phonons. 
Our neutron scattering data clearly reveal LA-|-TO^LO 
scattering processes, and an anomalous TO branch at 
F displaying a "waterfall" effect, as well as a strongly 
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temperature-dependent double-peak structure. 

A single crystal of PbTe {m ~ 25 g) was grown by 
a modified Bridgman technique, and characterized with 
x-ray and transport measurements. Details arc given 
in [16]. Hall measurements on the crystal yielded a 
temperature independent carrier concentration of 1.7 x 
10"'^'' electrons cm~^ from 4K to 300 K, and a room tem- 
perature mobility of 1540cm^/Vs, attesting to the high 
quality of the crystal. Inelastic neutron scattering mea- 
surements were performed with the time-of-flight Cold 
Neutron Chopper Spectrometer (CNCS) at the Spalla- 
tion Neutron Source, and (on the same sample) with the 
HB3 triple-axis spectrometer at the High Flux Isotope 
Reactor, both at Oak Ridge National Laboratory. Details 
are given in [16]. In the CNCS measurements, multiple 
datasets were acquired for different orientations of the 
crystal, and were combined in software to map the four- 
dimensional scattering function, S{Q,E), as a function 
of energy transfer E and wavevector transfer Q = q -f- r, 
with q a phonon wavevector and r a reciprocal lattice 
vector. Our data span multiple Brillouin zones, with 
larger coverages for higher incident neutron energies, Ei. 
Phonon dispersions were computed from density func- 
tional theory (DFT) using the linear response approach, 
using the same methods as in [17]. Using the phonon 
energies and polarization vectors from the DFT compu- 
tations, we calculated the dynamical structure factor for 
the scattering of neutrons by phonons. Details are given 
in [16]. 

One of the striking results of our investigation is the 
observation of a previously unreported "avoided cross- 
ing" between LA and TO phonon dispersions, associated 
with an anharmonic repulsion between these modes. In 
contrast with the harmonic dispersions calculated with 
DFT, the measured LA and TO branches repel each other 
strongly around q — (0,0,1/3). Fig. 1-a is for q along 
[OOL] in the {H,K,L) = (113) Brillouin zone at 300 K 
(where (H,K,L) denote reciprocal lattice units, rlu). The 
white lines plotted on top of the INS data are the har- 
monic phonon dispersions from DFT. Note that, accord- 
ing to our harmonic DFT calculations, the LA and TO 
branches should cross at q = 0.4 (we define q = L — 3 
in this zone) and E = 7.5 meV (in agreement with cal- 
culations of [18]). Measured LA phonon dispersions are 
linear for g ~ 0, with E'la increasing from at F to a 
maximum ^ 8mcV at q = 0.5, then curving back down 
to ~ 4.5 meV at the zone boundary (X), in good agree- 
ment with the DFT curves (see data in Fig. 1, panels 
a,c, and corresponding schematics labeling branches in 
panels b,d). However, the behavior of the TO branch 
is quite different from the DFT calculations in the har- 
monic approximation. The transverse optic (TO) branch 
has a minimum Eto — 4meV at F, and a steep disper- 
sion to ~ lOmcV at q ~ 0.2 rlu, and does not cross the 
LA branch. Instead, we observe a strong repulsion be- 
tween the two branches. The TO branch quickly rises 
above the top of the LA branch, while the intensity of 
the LA branch is strongly suppressed around q = 1/3. 



We measured the same effect in all other zones where 
the TO mode was observable ((H,K,L) odd). We note 
that the LA-TO crossing would normally occur between 
non-interacting modes of different symmetries along the 
symmetry lines, and thus the observed avoided crossing 
must arise from an anharmonic interaction. 

PbTe is an incipient ferroelectric material [19, 20]. 
While paraelectric when undopcd, it becomes ferroelec- 
tric upon substitution of a few percent of Pb by Sn or 
Ge [20]. The ferroelectric transition is of the displacive 
kind. It arises upon cooling through a critical temper- 
ature (Curie temperature, Tc), at which the "soft" TO 
phonon at the Brillouin zone center (F) slows to zero fre- 
quency, freezing in a lower-symmetry ferroelectric struc- 
ture [20, 21]. The TO phonon branch in PbTe behaves 
as in ferroelectrics, bending to low energies at F, and in- 
creasing in energy with increasing T [13, 14]. However, 
the TO mode never becomes fully soft, and the ferroelec- 
tric distortion is avoided [14, 20]. Previous first-principles 
computations of the electronic structure and phonons re- 
ported a strong sensitivity of the TO frequency on the 
volume of the system [17, 18, 22], and noted a possible 
coupling with the longitudinal acoustic (LA) phonons, 
which are compression waves [17]. Our measurements 
directly confirm this prediction, and show that the inter- 
action is strong at 300 K. 

The LA-TO interaction seen here is reminiscent of 
the TA-TO interaction in ferroelectric materials, such 
as perovskite PbTiOa [23]. However, the TO interac- 
tion with LA rather than TA modes is unique to PbTe, 
to our knowledge, and implies a strong anharmonic mix- 
ing, which will not occur at the harmonic level. The 
repulsion has its origin in an anharmonic coupling term, 
14,nhj which renormalizes the harmonic phonon frequen- 
cies most strongly at q where bare harmonic dispersions 
would otherwise cross, and is expected to have the form 
X I(LA|T4nh|T0)|V(SLA -Eto)- This mixing also 
affects the INS intensities, resulting in the LA extinction 
at q= 1/3, as we show using mode coupling theory [16]. 
The effect is observed even more strongly with increasing 
temperature (see data for T = 600 K in [16]), with the 
TO branch raising to higher E, and with an additional 
region of "collapsed" intensity at the bottom of the LA 
branch, close to X. The increase in the LA-TO interaction 
with T is expected, since it results from an anharmonic 
coupling, and the interaction will be stronger as phonon 
displacement amplitudes increase with T. 

Figure 1-c shows an interaction between LA, LO, and 
TO branches at the top of the LA branch (in the (002) 
zone). The LA dispersion is in good agreement with the 
DFT results and. in contrast with the (113) zone, its in- 
tensity is not extinguished at g = 1/3 (we define q = L—2 
in this zone). The absence of LA extinction can be re- 
lated to the vanishing TO mode structure factor in (002) 
(see [16]). The LO branch is seen clearly for 0.5 < q < 1, 
but is weak for < < 0.5, again in agreement with the 
dynamical structure factor. At q = 0.5 (both L = 2.5 and 
the symmetry-equivalent L = 3.5), some extra intensity 
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FIG. 1: CNCS data for PbTe at 300 K, showing the avoided-crossing behavior of LA and TO phonon branches in (a), the 
LA+TO— >LO scattering in (c), and the "waterfall" effect for the TO branch at F in (e). Solid and dashed white lines in 
(a,c,e) are harmonic dispersions calculated with DFT. (a,c): S{Q,E) data for dispersion along [0,0, L] in (113) and (002) zones, 
respectively. Red arrows in (a) point to the region of the avoided crossing of LA and TO modes. Red boxes in (c) show extra 
scattering intensity between LA and LO branches, (e): dispersion along [H,H,0] in (113) zone, showing the "waterfall" at F. 
(b,d,f) are schematics of the dispersions (blue lines), with blue rectangles representing diffuse extra scattering, and the bare 
TO branch as a thin black line. In all panels, pink diamonds indicate the position of the peaks in the TO scans at F (113), 
obtained with IIB3. Yellow square in (a,e) is Eto = 3.9 ± 0.2 meV from [13]. Data in (a,e) were collected with Ei = 12.1 meV 
(phonon creations). Data in (c) were collected with Ei — 3.5 meV (phonon annihilations, corrected for detailed balance). 
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is clearly seen between the LA and LO branches (it was 
also observed in other even zones), as indicated by red 
boxes in Fig. 1-c. We associate it with a 3-phonon inter- 
action process. The energy difference between LA and 
LO branches at g = 0.5 is 4meV, the same as the energy 
of the TO mode at F. Thus, we can write: LA(q = 0.5) 
+ TO{q = 0) O LO(q = 0.5), while conserving phonon 
energy and momentum [24]. The anharmonic LA-TO 
interaction around q — 1/3 is also reflected in the LA 
branch in (002), as a softened and broadened dispersion 
at low q. A more detailed analysis of the LA branch in 
this region is presented below. 

Fig. 1-e shows 5(Q, E) along [HHO] measured in the 
(113) Brillouin zone, and reveals a "waterfall effect" in 
the TO branch at F, also arising from the LA-TO interac- 
tion. The LA branch is very weak in this zone, but the TA 
and TO branches are clearly observed, in agreement with 
calculations [16]. The LA becomes visible only close to 
the X point (1,1,2), where it merges with the TA. While 
the TA branch in our measurement is in good agreement 
with previous reports and calculations [13, 17, 18], the 
TO branch is anomalous at the zone center. The branch 
gradually decreases to £' ~ 6.5 meV with q getting close 
to F, but at F the dispersion abruptly falls to very low E. 
We could not determine an actual lower bound, but in- 
tensity was detected down to ~ 2.5 meV, where it merges 
with the acoustic branches. This behavior is similar to 
the "waterfall" effect observed in ferroelectric materials, 
where the strong TA-TO interaction leads to an anoma- 
lous TO dispersion near F [25, 26]. In zones where the 
acoustic branch has a vanishing intensity, the anomaly 
is then expected to occur only at F [25], as we observe. 
This behavior was seen in both (111) and (113) zones. 
As a result, the scattering intensity at F shows a very 
broad distribution with E. We point out that according 
to our results, the TO branch dips down to low energies 
as Q F as generally expected for an incipient ferro- 
electric, but the LA-TO interaction introduces a complex 
profile of scattering intensity at F, extending to very low 
E. 

Fig. 2 shows data for constant-i? cuts along [0,0, L] in 
the (113) zone from measurements on both HB3 and 
CNCS. The two measurements are in excellent agree- 
ment, and the avoided crossing of LA and TO branches 
may clearly be seen. These data confirm our results 
shown in Fig. 1-a/c. We stress that the anharmonic LA- 
TO coupling is quite extended in reciprocal space. It 
is responsible for the waterfall effect at q = (in odd 
zones), the extinction of LA modes and avoided crossing 
with the TO branch around q — 1/3, and the diffuse scat- 
tering intensity at the top of the LA branch {q — 0.5). 
This large phase-space for the coupling, together with the 
strong coupling, is expected to affect the contribution of 
LA modes to the thermal conductivity of PbTe. Since 
LA modes are fast-propagating phonons (large group ve- 
locities for q < 0.5), they are expected to be important 
heat carriers, and thus the LA-TO coupling reported here 
probably plays an important role in explaining the low k 
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FIG. 2: Constant-i? intensity cuts along [0,0, L] in the (113) 
zone, showing the avoided crossing behavior of the LA and TO 
branches. Data were obtained with two different instruments: 
HB3 (a) and CNCS (b). Both datasets are for the same 
PbTe crystal at 300 K. Stars (resp. dots) serve to identity 
the TO (resp. LA) mode positions. The TO mode (stars) is 
always "inside" the LA mode, indicating that the two phonon 
branches do not cross. 



of PbTe. 

We investigated the LA branch in the (002) zone in 
more details, by performing constant-Q cuts of the CNCS 
data, and fitting the position of the LA peak (Gaussian 
fits, corrected for instrument resolution). The results are 
shown in Fig. 3. The LA dispersion shows an anomalous 
dip around q ~ 0.2, and the phonon linewidth F exhibits 
a peak at the same wavevector, reaching above 2 meV 
at the maximum. Scans measured on HB3 at (0,0,3.8), 
equivalent to q = 0.2, also showed a very broad distribu- 
tion, while clear LA peaks could be observed at q = 0.1 
and q = 0.3. These anomalies can be directly related to 
the anharmonic interaction with the TO branch, which 
results in a down-shift and an increased damping of the 
LA phonons. We note that in the range q ~ 0.15 — 0.35, 
the shape of the LA peak is broad and complex, with ex- 
tra broad scattering intensity detected between the LA 
and TA peaks (this was observed in both our CNCS and 
HB3 measurements). The data for the cuts are shown in 
[16]. 

Investigations of the INS intensity at F reveal a broad 
and complex energy profile for the TO mode, as well as 
a strong dependence on temperature, indicative of an- 
harmonicity. Scans at F =(1,1,1) and (1,1,3) measured 
with BBS at multiple temperatures are shown in Fig. 4. 
The TO intensity is broad in E, in agreement with the 
CNCS data, and it is also strongly dependent on T. At 
T = 100 K, a single well-defined peak is observed at 
E — 3.2 meV, with a broad tail at higher E. This peak 
position agrees well with i?TO reported in [14]. With in- 
creasing T, the peak shifts to higher energy, as occurs in 
ferroelectrics, and in agreement with previous reports in 
[14]. Besides this shift, we observe that the INS intensity 
profile becomes more complex, with a second peak ap- 
pearing at i? ~ 6 meV, while the low-E peak is strongly 
broadened, corresponding to damping of the TO mode. 
A double-peak structure is observed at 300K, especially 
clearly in HB3 data at (111) - and also in CNCS data 
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FIG. 3: LA mode phonon dispersion (-Ela) and linewidth 
(Fla, corrected for instrument resolution), showing anomaly 
at the phonon wavevector q ~ 0.2 along [0,0,L]. The data 
are results of fits of constant-Q cuts through CNCS measure- 
ments with Ei — 3.5 meV (phonon annihilation processes). 
The energy resolution (FWHM) was 0.05 meV at i5 = meV, 
increasing to 0.2 meV at i5 = 8meV. The q resolution was 
about 0.02 rlu. The red dashed line is a guide for the eye 
(E — 7.9 X sm{nq)). Error bars for Ela are comparable to 
size of markers. 



at (113)-, and is consistent with the reflectivity data of 
[27] (reproduced here for comparison). The high-i? peak 
becomes more intense, and also shifts up in energy with 
increasing T. We fitted the positions of the two peaks 
with Gaussians, and we find that the low-E peak corre- 
sponds to the calculated bare i?TO at F, while the high-i? 
peak corresponds roughly to the bottom of the renormal- 
ized branch, as seen in Fig. 1 (pink diamonds) . More de- 
tails are given in [16]. We note that since the TO mode 
is sensitive to strains, different local environments could 
lead to a broad distribution of TO frequencies. Thus, it 
is possible that fluctuations (such as reported in [28]), 
or anharmonic displacements of large amplitudes, could 
produce the broad TO spectrum. 

Our inelastic neutron scattering measurements reveal 
a strong, unprecedented coupling between LA and TO 
phonon branches, in particular along the [001] direction. 
This coupling results in an avoided crossing behavior for 
the measured dispersions around q ~ 1/3, with a soft- 
ening and damping of the LA branch and a repulsion 
of the TO branch. We also observe some intensity in 
S{Q,E) that can be associated with three phonon pro- 
cess LA-hTOoLO at [0, 0,q = 0.5] (involving the TO 
mode at F), and a waterfall behavior for the TO branch 
at r{q — 0). These observations are all consistent with 
an extended anharmonic interaction between the LA and 
TO phonons in PbTe, in agreement with DFT calcula- 
tions [17]. The large phase-space for the coupling, to- 
gether with the coupling strength, are thus expected to 
significantly affect the contribution of LA modes to the 
thermal conductivity of PbTe. The effects described here 
for PbTe may also be relevant to explain why materials 
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FIG. 4: Profile of the TO mode in energy at the zone center, 
measured with HB3 and CNCS, showing the broad double- 
peak structure and its change with temperature. The neutron 
scattering intensity (squares and diamonds) vs E were mea- 
sured in constant-Q mode at F points (113) in (a) and (113) 
in (b), as a function of temperature. Data for different tem- 
peratures are offset vertically for presentation. Continuous 
lines are fits of the data with two Gaussians for the double- 
peak structure, and an additional Gaussian for the elastic line. 
Panel (b) also shows the reflectivity data of Burkhard et al, 
for comparison [27]. 



crystallizing in the sodium chloride (rocksalt) structure 
typically have lower thermal conductivities than materi- 
als with cesium chloride structures [29]. We also point 
out that the LA and TO modes along [100] involve vi- 
brations of Pb and Te atoms with displacements parallel 
to [100] directions, and are thus possibly related to the 
fluctuations recently reported in [28] . The TO mode (fer- 
roelectric mode) interacts with acoustic phonons over a 
wide range of frequencies, and bears some resemblance 
to a rattling mode, which is a consequence in this case 
of a strong anharmonic coupling. This provides a mech- 
anism for a material to behave as if it has a very soft 
lattice for the purpose of heat conduction, while at the 
same time not being unduly soft from the point of view of 
stiffness or melting, thus providing an avenue for finding 
very low lattice thermal conductivity crystalline materi- 
als that nonetheless are stable at high temperatures. The 
present results provide a microscopic picture for the low 
thermal conductivity of PbTe and point towards a new 
direction for finding good thermoelectrics. 
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A. Sample Preparation and Characterization 

A stoichiometric mixture of high purity Pb (99.999%) 
and Te(99.9999%) shot (20.72 g Pb, 12.76 g Te) was 
loaded into a carbonized, round-bottom silica tube, that 
was then evacuated and sealed. The mixture was heated 
to 1050° C at 2°C/min, homogenized and mechanically 
mixed for 30 h at temperature, and then cooled to room 
temperature at 5°C/min. The resulting boule was re- 
moved from the silica ampoule and any slight oxidation 
lightly sanded from the surface. The boule was broken 
into mm size pieces, transferred to a silica Bridgman cru- 
cible, evacuated and sealed. The Bridgman ampoule was 
placed in a box furnace with a natural vertical tempera- 
ture gradient, heated to 1050°C at 2°C/min, soaked for 
24 h, cooled over several hours to 980° C, slowly cooled 
(0.6 - l°C/h) to 750°C, followed by a l°C/min cool to 
room temperature. The resulting single crystal was re- 
moved from the silica ampoule. Small slices were cut 
from the crystal using a slow speed diamond saw for 
additional x-ray, chemical and transport characteriza- 
tion. Hall measurements on the crystals yielded a tem- 
perature independent carrier concentration of 1.7 x 10^'' 
electrons/cm'^ from 4K to 300 K, and a room tempera- 
ture mobility of 1540cm^/V s. (The mobility at 4 K is 
an astounding 1.7 x 10^ cm^/ V s). 

B. Inelastic Neutron Scattering Measurements 

In the measurements on CNCS, the sample was 
mounted with an (HHL) plane approximately horizon- 
tal, at r = 300, 500, 600 K using a low-background fur- 
nace. In CNCS measurements, multiple datasets were 
acquired for different orientations of the crystal, corre- 
sponding to rotations about the vertical (//[— 110]) axis 
(typically in steps of 2 or 4°). The datasets for differ- 
ent rotations were then combined in software to produce 
a four-dimensional sampling of the scattering function 
iS'(Q,i?) of the crystal. Small departures from the exact 
horizontal (HHL) orientation were corrected in software 
when analyzing the data. Several incident energies were 
used, Ei = 3.51, 12.03, 25 meV, with corresponding en- 
ergy resolutions (full width at half maximum) at the elas- 
tic line 2aE = 0.05, 0.5, 1.0 meV. At each incident energy 
and each temperature, multiple datasets were acquired 
for different orientations of the crystal, corresponding to 
rotations about the vertical axis (typically in steps of 2 
or 4°). The datasets for different rotations were then 
combined with the software MSLICE to produce a four- 
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FIG. 1: CNCS S{Q,E) data for dispersion along [0,0, L] in 
(113) Brillouin zone, showing the avoided-crossing behavior 
of LA and TO phonon branches at 300 K in (a), and at 600 K 
in (b). White lines are the harmonic phonon dispersions cal- 
culated with DFT. The pink diamonds indicate the positions 
of the peaks in the TO scans at F (113), obtained with HB3 
at 300 K and 670 K. 



dimensional sampling of the scattering function S'(Q, E) 
of the crystal [1]. These data were then shced to pro- 
duce scattering intensity as a function of energy transfer 
E and momentum Q parallel to a chosen direction. 

Figure 1 shows the CNCS S{Q, E) data for dispersion 
along [0,0,L] in (113) at 300 K and 600 K. As may be 
seen on the figure, the avoided crossing of LA and TO 
branches (see main text) is observed more strongly at 
T = 600 K than at 300 K, with the TO branch raising to 
higher £', and with an additional region of "collapsed" 
intensity at the bottom of the LA branch, close to X. At 
600 K, the TO branch shows only one maximum of inten- 
sity at the F point, in agreement with the HB3 constant- 
Q scans. 

In order to investigate the longitudinal acoustic (LA) 
phonon dispersion in the (002) zone. We took some 
constant-Q cuts from the CNCS data in Fig. 1 (main 
article), integrating of a range of AQ[ooi] = 0.03 rlu at 
each step (AQjug] = 0.15, and AQj.hq] — 0-15 rlu). 
Resulting data are shown in Fig. 2. The data in the 
figure were corrected for detailed balance, and are plot- 
ted vs positive phonon energy. The LA peak is clearly 
observed, and one may notice that it is sharp at low-q, 
broadens around q = 0.2, and then sharpens again to 
the top of the branch at q = 0.5. The LA peak posi- 
tion and peak width was fitted with a Gaussian, and the 
phonon linewidth was obtained as the peak full-width 
at half maximum, after correction for instrument resolu- 
tion (including the effect of Q-resolution and dispersion 
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FIG. 2: (a) Constant-Q cuts for q along [OOL] in (002) zone 
from CNCS data at 300 K {Ei = 3.5 meV). Curves are ver- 
tically offset for clarity. Dashed lines are guides for the eye, 
indicating the positions of the transverse acoustic (TA) and 
longitudinal acoustic (LA) phonon peaks. The anomaly in the 
shape and width of the LA peak is clearly seen in the range 
g = 0.15- 0.25 rlu. 

slope on the E resolution). Results are plotted in the 
main manuscript. We note that the linewidths were dif- 
ficult to extract below q = 0.05 (because of nearby TA 
branch) and above q = 0.45 (because of very broad LA- 
LO interaction), and are not reported. The transverse 
acoustic peak (TA) is also observed, although weakly (in 
agreement with dynamical structure factor calculations, 
see below), and its position is in agreement with results 
published in [2]. In the region of the LA-TO avoided 
crossing [q — 0.1 — 0.3), some broad scattering intensity 
is observed between the LA and TA peaks, which may 
also be associated with the strong anharmonic behavior 
in this region. 

Additional INS measurements were performed with 
the triple-axis neutron spectrometer HB-3 at the High- 
Flux Isotope Reactor at Oak Ridge National Labora- 
tory. We used pyrolitic graphite monochromator and 
analyzer (in (002) reflection condition), and worked in 
constant final energy mode, with Ef — 14.7 meV. The 
horizontal coUimation angles were 48, 40, 40, and 120 
minutes of arc for the SoUer collimators prior to the 
monochromator, between monochromator and sample, 
between samples and analyzer, and between analyzer 
and detector, respectively. An (HHL) crystallographic 
plane of the sample was aligned with the horizontal scat- 
tering plane. The sample was mounted either inside 
an helium cryostat or a furnace for measurements at 
T = 10, 50, 100, 150, 200, 250, BOOK (cryostat), or T = 



300, 500, 670 K (furnace). The signal from the empty fur- 
nace was also measured. A vanadium rod was measured 
to determine the energy resolution of the instrument at 
the elastic line (zero energy transfer). The profile was 
Gaussian, centered on E = OmeV, with a l.lmeV full 
width at half maximum. 

Scans at F =(1,1,1) and (1,1,3) were measured with 
HB3 at multiple temperatures (see main article). The 
IIB3 data were fitted with Gaussians for the low- and 
high-i? peaks (and another Gaussian for the elastic peak). 
The positions of the centers of the two Gaussians are plot- 
ted as a function of temperature in Fig. 4. The widths 
of the excitations were extracted by corrected for the in- 
strument resolution function, which was calculated as a 
function of energy transfer, and for the instrument condi- 
tions used. The peaks observed in I{E) are significantly 
broader than the instrument resolution (see data in main 
article) . 

We plot our results for the positions of the two peaks 
as a function of temperature in Fig. 3-a. As can be seen 
on this figure, the position of the \ow-E peak at 300 K 
coincides with the value of -Eto a-t 296 K previously re- 
ported by Cochran [2] . The value reported in [3] for -Eto 
at 295 K is slightly smaller than our measurement (by 
0.4 meV), but still in good agreement, considering that 
-Eto a-t the zone center could be sensitive to doping lev- 
els. The T-dependence of Eto,i is comparable to that 
reported in [3], although we find energies that are con- 
sistently slightly higher. The energy of the second TO 
peak, Exo,2: is about 2meV higher than Exo,!; at all 
temperatures. Both Eto.i and Eto,2 follow a behav- 
ior reminiscent of that in ferroelectrics, but we could not 
successfully fit their T-dependence to the typical ferro- 
electric mode behavior (E^ (x {T ~ Tc)) over the full 
temperature range. Also, we observe an apparent sat- 
uration in the behavior of Eto,i and Eto,2 at high T, 
which could be related to the saturation in amplitude of 
the structural flutuations reported in [4] . 

Our results for the widths of the peaks are shown in 
Fig. 3-b. A significant increase of the energy widths of the 
TO peaks (FWHM with of Gaussian peak fits, corrected 
for resolution), F, is observed with increasing tempera- 
ture for both TO peaks. The increase in F is related 
to a damping of phonon excitations, and is generally ex- 
pected for vibrations in anharmonic potentials. As may 
be seen in Fig. 3, the width (FWHM) of the low-energy 
TO peak (Fxo.i) is more strongly dependent than that 
of high-energy TO peak (rTo,2)- Between 50 and 500 K, 
Fto.i nearly doubles. At 670 K, we could not extract 
the position or width of this peak anymore, as it became 
too broad and the Exo,2 became more intense (see main 
article) . The precise T-dependence of phonon linewidths 
in solids is a topic of current interest, and definite trends 
remain to be established [5]. Both Fto,i and Fto,2 ap- 
pear to saturate (or at least increase more slowly) as 
T reaches above 300 K, which is rather unexpected, but 
could be related to the saturation in the amplitude of 
fluctuations reported in [4]. It is also interesting to com- 
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FIG. 3: Position (a) and linewidth (b) of peaks fitted to HB3 
data at the (111) and (113) zone centers. 



pare the position of the two TO peaks at F with the 
dispersion branches reported in [2] and the INS intensity 
measured with CNCS (main article). The position of the 
peaks observed at 300 K in HB3 data is shown with pink 
squares in Fig. 1 (superimposed on the S{Q, E) measured 
with CNCS). We can see that i?TO,2 matches the top of 
the region of diffuse scattering at the zone center, and 
corresponds to the energy where the dispersion would be 
expected to converge for Q — > 0. On the other hand, 
-E'TO,! matches the position of TO mode in the harmonic 
DFT calculations. 



C. Dynamical Structure Factor and Mode 
Coupling Calculations 

The phonon frequencies and the polarization vectors 
were calculated on a fine 40x40x40 grid of fc-points in 
reciprocal space using the linear-response approach, fol- 
lowing methods in [6] . The lattice parameter used in the 
calculations was a = 6.376A, which gave dispersions in 
close agreement with the measured dispersions at 30GK, 
in particular for the energy of the TO branch close to F 
[2, 6]. From the calculated phonon frequencies and po- 
larization vectors, we computed the dynamical structure 
factor in [110] and [001] directions in the (113), (004), and 
(220) Brillouin zones. The neutron scattering intensity 
was calculated based on the DFT phonon polarizations 
Cd^s and energies Es , where d is an index for the atom in 
the crystal basis (Pb or Te), and s is the phonon mode 



index, for phonon wavector q: 

S{Q,E) cx ^^-l|^^exp(*Q.r,)exp(-W'rf) 



(Q • I X {n, + l)d{E - Es)6{Q - q - r) , 

with r a reciprocal lattice vector, Q the wavevector 
transfer, the position vector for atom d inside the 
unit cell, AI^ its mass, and its coherent cross-section 
[7]. The thermal occupation function is taken into ac- 
count through the Bose-Einstein distribution Ug (we used 
T = 300 K in our calculations). The Debye- Waller fac- 
tor Wd was also calculated from the polarizations and 
energies: 



Wd = 



4Md 



E 



|Q • ed,. 

e7 



Results are shown in Fig. 4. The calculated dynamical 
structure factor is in good agreement with the measure- 
ments along [001] in the (002)- (004) zones, with only the 
LA and LO branches having intensity, and the inten- 
sity of the LO becoming stronger as it gets closer to the 
zone boundary. The main difference is the extra scatter- 
ing intensity observed experimentally halfway between F 
and X between LA and LO, as discussed in the main 
text. Similarly, along [110] in the (113) zone, only the 
TA and TO branches have strong intensity, as observed 
experimentally, although the the harmonic calculations 
do not capture the waterfall effect in the TO branch at 
F. The structure factor calculation for harmonic, non- 
interacting phonons show the strongest discrepancy with 
measurements for the [001] dispersion in the (113) zone, 
as discussed in the main text. Although the TO, LA, and 
TA are all visible, the experiment shows a gap between 
the LA and TO (avoided crossing) around g = 1/3, where 
harmonic calculations predict overlapping intensity. This 
discrepancy is associated with the LA-TO interaction, 
and we perform a simple model calculation below to show 
how this interaction can lead to the extinction of the LA 
branch near the avoided crossing. 

We performed simple mode coupling calculations, to 
investigate the extinction of the LA branch observed in 
our data along [001] in the (113) zone. We limited our- 
selves to a simple one-dimensional model with only LA 
and TO branches, and followed the model described in 
[8]. The non-interacting harmonic dispersions were ob- 
tained from the DFT calculations described above. We 
investigate the effect of a simple g-independent coupling, 
A, on the dispersions and the scattering intensity. The 
coupling introduces non-diagonal elements in the dynam- 
ical matrix: 



A* 



A 

2 

TO 



where WLA(g) and ujTo{q) are the bare frequencies of 
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the LA and TO branches at wavcvector q, obtained with 
DFT. We neglect non-diagonal terms in the damping ma- 
trix, 



Tla 
Tto 



where Fla and Fto are the linewidths of LA and TO 

modes, which wc take cqnal to 1 meV and 2 mcV, respec- 
tively. At high temperature T, the scattering intensity 
is calculated as I{q,uj) cx kT/oj f{q)* luiG{q,u]) f{q), 
with G{q,Lj) = {V{q) - iujr{q) - co"^ E2)~\ with E2 a 
2x2 unit matrix, and f{q) a 2- vector defining the rel- 
ative strengths of LA and TO branches, which we keep 



equal to [1,1] for simplicity. We show the results for 
A = OmeV^ and A = 9meV^ in Fig. 5. 



The interaction produces a repulsion between the 
branches, and an "avoided crossing" where the bare 
phonon branches would normally cross. The interac- 
tion also produces the weakening in intensity of the LA 
branch near the avoided crossing, as observed experimen- 
tally. We note that our simplistic model with a constant 
A(g) would tend to produce unstable LA phonons near 
r, but we cannot entirely trust such a simple calculation. 
In particular, the real interaction, besides being three- 
dimensional, could have a non-trivial q dependence. 
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FIG. 4: Dynamical structure factor (color maps) of PbTe for harmonic phonons, calculated from phonon polarization vectors 
obtained by linear-response DFT (see text). White curves are harmonic phonon dispersions from linear-response DFT calcu- 
lations, a: dispersion along [001] in (004) zone, b: dispersion along [110] in (004) zone, c: dispersion along [001] in (113) zone, 
d: dispersion along [110] in (113) zone, e: dispersion along [001] in (220) zone, f: dispersion along [110] in (220) zone. 
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FIG. 5: Calculation of effect of LA-TO mode coupling on dispersion branches and scattering intensity. (a,b): dispersions, 
without and with coupling, respectively. (c,d): scattering intensity, without and with coupling, respectively. In both (b) and 
(d), the strength of the coupling is A = 9meV^. 



